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Abstract 

Extending our previous work on 2D growth for the Laplace equa- 
tion we study here multidimensional growth for arbitrary elliptic equa- 
tions, describing inhomogeneous and anisotropic pattern formations 
processes. We find that these nonlinear processes are governed by an 
infinite number of conservation laws. Moreover, in many cases all dy- 
namics of the interface can be reduced to the linear time-dependence 
of only one "moment" Mq which corresponds to the changing volume 
while all higher moments, Mi, are constant in time. These moments 
have a purely geometrical nature, and thus carry information about 
the moving shape. These conserved quantities (eqs. (7) and (8) of this 
article) are interpreted as coefficients of the multipole expansion of 
the Newtonian potential created by the mass uniformly occupying the 
domain enclosing the moving interface. Thus the question of how to 
recover the moving shape using these conserved quantities is reduced 
to the classical inverse potential problem of reconstructing the shape 
of a body from its exterior gravitational potential. Our results also 
suggest the possibility of controlling a moving interface by appropriate 
varying the location and strength of sources and sinks. 
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PACS numbers: 47.15. Hg, 68.10.-m, 68.70.+W, 47.20. Hw. 
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Many seemingly different pattern formation processes have much in com- 
mon, both in their mathematical description and in their physical behav- 
ior. Among them are the famous Stefan problem (freezing of liquid), flow 
through porous media, the Rayleigh-Taylor instability, electrodeposition of 
metals, corrosion, combustion, growth of bacterial colonies, dynamics of 
earth cracks, diffusion-limited aggregation(DLA), etc. The common feature 
shared by these processes is the existence of an evolving interface. The prob- 
lem of the evolution of the interface in these processes attracts considerable 
attention (see, for example, the book ||]) both because of its great practi- 
cal importance and because of its connections with such fields as dynamical 
chaos, nonequilibrium physics, and fractal growth (see DLA 0). 

A general scheme for these processes is as follows. There is a linear PDE 
(frequently of the second order) for the scalar field determining the process. 
For example, this is the diffusion equation for the Stefan problem and the 
Laplace or Helmholtz equations for electrodeposition. This scalar field is 
temperature in the Stefan problem, pressure in flows through porous media, 
concentration of the nutrient for bacterial growth, electrostatic potential in 
electrodeposition, probability of the next jump in DLA, etc. Appropriate 
boundary conditions are imposed both on the moving part of the boundary 
(the interface) and non- moving part of the boundary (the outer walls). In 
addition, a law of interface motion is given in terms of the local behavior 
of the main scalar field. (Typically the local velocity of the interface is 
proportional to the gradient of the scalar field near the interface) . The main 
question is: 'what is the evolution of the interface?' 

It is remarkable that some of the problems mentioned above are ex- 
actly solvable in two dimensions ||-[||] despite the nonlinearity of these 
processes. These problems were solved with the help of time-dependent 
conformal mapping which cannot be extended to 3D except for a few trivial 
cases. In 2D, it was found that such processes as two-phase flows in porous 
media, electrodeposition, and slow solidification in a supercooled liquid or 
from a supersaturated binary solution are governed by an infinite number 
of the constants of motion which were obtained explicitly in several spe- 
cial cases ||— 1|. These constants of motion are related to the conserved 
moments proposed by Richardson ||. This invariance is quite subtle and 
disappears when realistic physical perturbations such as surface tension or 
random noise are added. 

In the 3D case very few exact analytical results are known: a constant- 
velocity paraboloid RJ, and a self-similar ly growing ellipsoid The only 
known way to obtain these solutions is by using one of the eleven coordinate 
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system in which the 3D Laplace equation is separable (9| . One then considers 
level surfaces as moving interfaces. It is clear that this method does not work 
when the shapes are time-dependent. The traditional attitude is that the 
main obstacle in obtaining nonperturbative exact results in 3D is the lack 
of (nontrivial) conformal mappings unlike in the 2D case. But is not this 
statement too strong? 

This article is a natural extension of previous work || to more general 
and realistic multidimensional growth problems; and not only for the Laplace 
equation as was done in |J-§|, but also for arbitrary elliptic equations de- 
scribing for example, inhomogeneous and anisotropic diffusion in solidifica- 
tion, inhomogeneous dielectric functions and screening in electrodeposition, 
and inhomogeneous viscosity for flows through porous media. It turns out 
that these nonlinear processes also possess remarkable properties (an infinite 
number of conservation laws) similar to the ones mentioned in |J, || for the 
2D Laplacian case; and these properties do not depend on the dimension 
of the process considered. Thus we show that, contrary to the traditional 
attitude, we do not need a conformal mapping for this invariance. Rather 
this invariance originates from the more general property: the elliptic nature 
of the equation for the scalar field. 

Let us state now the following D-dimensional problem: 

L(u) = div(p(r) grad-u(r)) + q(r)u(r) = (1) 

for r £ B C R D where the domain B is bounded by the nonmoving exterior 
boundary, E, and by the moving interior boundary T(t) (t is time) which 
is the interface separating the domains B and A. An interior domain A 
contains the origin and is surrounded by the moving interface T(t) as shown 
at the Fig.l. Here p and q are given functions of r = (x%,X2, ...,xr>)- The 
boundary conditions imposed on u are: 

cU| E = G(£) , (2) 

u(T(t)) = . (3) 

The left-hand side (LHS) of (2) means the normal component of grad(u) 
evaluated at S. 

There can also exist point-like sources and sinks in the domain, B, lo- 
cated at rfc and having strengths, Sk, (k = 1, 2, N), so that near r& , w(r) 
diverges and satisfies 
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u = Sk/\r — rk\ D 2 + smooth function (4) 

if D > 2, or 

u = Sk log|r — r&|+ smooth function 

if D = 2. 

The law of motion of is 

u n = -p(r) 9 n «| r ( t ) , (5) 

where the LHS is the normal component of the velocity of T(t). 

Equations (1) and (5) together with the boundary conditions (2)-(4) 
complete the mathematical description of the motion of T(t). If, for example, 
D = 3, p(r) =const, and q(r) = 0, this describes slow solidification or two- 
phase flow in porous media. 

In this paper, I show that, in spite of the complexity and nonlinearity 
of the processes described by eqs. (l)-(5), these processes are governed by 
an infinite number of conservation laws. Namely, if the outer boundary, 
S, is very far from the origin there is an infinite set of numbers C[ (I = 
D — 1, D, D + 1, oo) defined as 



Ci = dMi/dt = d/dt{ / ^ d D r) (6) 



( 7 il U ' 

B 

which are conserved during the evolution of the hyper surface^], T(t), and 
equal 

7T D / 2 N 

Cl = 2 T(D/2) Y, k=1 s kP( r k)M^k) ■ (7) 

Here, T(n) is the Gamma function (not the interface T(t)) and tpi is the 
arbitrary solution of L(u) = which decays at infinity no slower than r~ l 
and having singularity only at the origin^. If g(r) = we have one more 
conserved quantity: Co, which is the rate of the change of the volume of B 
when L is Laplacian, and which satisfies 

1 We do not consider here the passing of the interface through the singularities. 

2 The functions ipi as well as quantities C; and Mi are labeled in general by more than 
one number (which is I here). See for example the Eq.(ll) below. But for simplicity and 
without the loss of generality we drop all labels except the I almost everywhere. 
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C = dM /dt = d/dt{[ p(r)d l 

JB 

tt d / 2 n r 
= 2 f^/2)E fe=1 s ^(r fc ) + y s G(S)p(S)dS , (8) 

Here we took tpo = 1 which is a solution of L(u) = when q(r) = 0. 

We think that the knowledge of C\ 's together with the initial Mi 's defined 
in (6) (the latter s are uniquely determined by the initial shape of the inter- 
face ) could describe the whole moving shape in many of important physical 
cases. We consider eqs. (7) and (8) as the main result of this work. 

Formula (7) follows immediately from the following considerations: 
Since 

/ ^id D v- f ^d D r= f ^v n dTdt 

JB(t+dt) JB(t) Jr(t) 



we have 



d/dt([ ^d D r) = f ^v n dT 
JB(t) JT(t) 



?(t) ./n/) 
Further, because of (5), it equals 



= / {-p^id n u)dT 
JT(t) 



■(f) 

and, finally, in view of (3) this expression is 



p(u d n ipi - tpi d n u) dT 

T(t) 

= / {y/pud n (y/pipl) - y/pipld n (y/pu))dr . (9) 

Jr(t) 

Applying Green's theorem to the functions ^Jpu and y/pipi we find that 
the RHS of eq. (9) is given by 

/ div( A /pugrad(- v /p , 0i) ~ y/P^l g^ad(y/pu)) d D r 

JB 

+ Jp(ud n ipi-ipid n u)cE (10) 

Ylk=i J p<yU dn ^ ~ ^ dn u ^ dTi • 
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Here the summation is over the point-like charges, Sk, and 7^ denotes the 
surface of the infinitesimal hypersphere around the r^. 
Considering the RHS of (10) one can see that: 

(i) the volume integral over B vanishes because of (1); 

(ii) the surface integral over £ also vanishes if the outer boundary, X, 
is far removed from the center and if ipi decays at infinity stronger than 
l/r D ~ 2 . (When ipi = 1, this integral is not zero but equals J s G(S)p(S) dT, 
as it is in the RHS of (8)); 

(iii) the contribution of the first integrand to the surface integral over 
the 7^ is zero, while the integral from the second term equals 

TT D ' 2 

since p(r) and ipiir) are regular near the r k and due to Gauss's theorem. 

Thus the RHS of (10) equals the RHS of (7) (or (8) when q(r) = 
and / = 0), so we have obtained the infinite set of the conserved quantities, 
Ci, if the sources and sinks are nonmoving (i.e. when s k and are time- 
independent). Moreover, if q{r) = and no point-like singularities are in 
the domain B (i.e. if all = 0), then all the dynamics of the interface T(t) 
has been reduced to the linear time-dependence of only one "moment", Mo, 
which is the volume of the phase B if L is Laplacian. All higher moments, 
Mi, are constant in time. Note also that all moments Mi have a purely 
geometrical nature, and thus carry information about the moving shape. 

It is worth mentioning the physical interpretation of the derived invari- 
ants, Mi, in the special case when the operator L is Laplacian: p = 1, q = OR. 
For D = 2, if the ipi are chosen as ipi = z~ l where z = x + iy, these integrals 
coincide with those previosly found in explicit form Q via the coefficients 
of the appropriate conformal map. These are analogs of the Richardson 
moments, whose invariance was found earlier for the interior Hele-Shaw 
problem || . 

In the 3D Laplacian case, one can choose the ipi to be a set of spherical 
functions: 

Here, r, 9, and <f> are the polar coordinates and (9) are the asso- 
ciated Legendre polynomials. In this case, the moments, M^ , are the 

3 In this case it was noted by Howison @ that the integral of a harmonic function over 
the infinite domain enclosed the moving bubble is a constant in time. 
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coefficients of the multipole expansion of the Newtonian potential at an ar- 
bitrary point of the empty interior domain A, if the potential is created by 
the mass uniformly occupying the domain B. Thus the question of whether 
one can recover the moving shape using only the numbers, Mi, introduced 
in (7) is now reduced to the classical inverse potential problem [ [l0| ] for the 
reconstruction of the shape of a body of constant density from its Newtonian 
potential. Our case corresponds to the exterior problem (where the poten- 
tial is given in the empty hollow of the body: in the phase A). The author 
believes that the connection between pattern formation studies and the in- 
verse potential problem is especially important and deserves close attention. 
We do not discuss this problem here. Rather we merely note that in 3D 
(unlike the 2D case) there is no description of a body (with the exception 
of the sphere) in terms of a finite number of nonzero moments, Mj. (For a 
detailed description of these difficulties see [11]^|). 



Note that it is also possible in the general elliptic case (when L is not 
Laplacian) to preserve the interpretation of the conserved quantities Mi as 
coefficients of the orthogonal expansion by choosing the Green's function 
G(r, ri) of the operator, L, to be the integrand in (7), since G(r, r{) satisfies 
the conditions imposed on the ipi if r G B and ri E A. Using the Green's 
function expansion and by choosing the orthogonal set of the eigenfunctions 
ipi of L(ipi) = bounded at infinity and divergent at the origin and the com- 
plementary set of eigenfunctions ipi of the same equation but with opposite 
asymptotics, we have 

U{r) = V j>i(v)d/dt(f ^ l {r 1 )d D r 1 ) . 

We remark that our main result (7) holds also for time-dependent 
and rfc (since we never used the time-independence of sj~ and in obtaining 
(7)). Although the Ci are no longer conserved, the problem is still integrable, 
since Ci are known functions of time (if the time-dependence of Sk and 
is given). Thus the moments, Mi, are easily controlled parameters, namely 
they are just primitives of the time-dependent RHS of (7). In this way, one 
might be able to govern the motion of the interface by the proper choice of 
the Sfc and r^. 

It should also be mentioned that there could be a few exceptions among 
the C;'s (only one in the 3D Laplacian case) for which eq. (7) is not valid and 
thus they may not be conserved. These nonconserved Ci's correspond to the 

4 I am grateful to M. Brodsky who brought my attention to this book. 
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V>z's decaying at infinity but not stronger than r^ 2 ~ D \ For the 3D Laplacian 
case the only nonconserved quantity is C\, which corresponds to tpi = r _1 . 
Since we do not know the time-dependence of C\ one could suppose that 
the description of the interface is now less complete. However, although C\ 
is not conserved, we have the conserved quantity Cq (see eq. (8)). In other 
words we think that the set {Mi, M2, M3, ...} describes the shape with the 
same degree of completeness as the set {M , M 2 , M 3 , ...}. This question 
concerning nonconserved quantities among the C/'s does not arise for the 
interior problem when the scalar field u is given in the phase enclosed by 
the moving interface T(t) (in our case in A instead of B). 

In conclusion we pose several questions that arise from these studies and 
which we beleive merit attention: 

1) Does the relation expressed by (7) really mean complete integrability 
of the multidimensional growth? 

2) If yes, to which nonlinear evolutionary PDE's do these constants C\ 
correspond? 

3) Is there a Hamiltonian structure for these systems? 

4) Do finite-time singularities (cusps?) exist here, as in the 2D case? 
And, if yes, how may surface tension regularise them? 

5) How can one recover the moving shape from the given set of moments, 
Mp. 
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